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A  Globally  Stable  Adaptive  Controller 
for  Multivariable  Systems 


R.  P.  Singh  and  K.  S.  Narendra 
Center  for  System  Science,  Tale  University 


Introduction 

Soon  after  the  problem  of  adaptively  controlling  a  single- input  single-output 
(SISO)  system  in  a  stable  fashion  was  resolved  In  1979  (1-3],  interest  shifted  to 
related  theoretical  problem.  One  of  the  principal  questions  currently  being 
investigated  is  the  possibility  of  extending  these  results  to  multivariable 
system  [4-10] . 

The  adaptive  control  problem  can  be  broadly  divided  into  two  parts  -  an 


algebraic  part  dealing  with  a  specific  parametrlzation  of  the  plant  and  an 
analytic  part  dealing  with  the  adaptive  laws  and  the  resulting  problem  of  con¬ 
vergence.  For  stable  adaptive  control  of  SISO  plants,  certain  assumptions  re¬ 
garding  the  plant  transfer  function  have  to  be  made.  In  particular  it  is 
assumed  that  the  designer  has  the  knowledge  of 

(i)  the  relative  degree  n*  of  the  plant  transfer  function, 

(il)  the  sign  of  the  high  frequency  gain  k  » 

P 

and  (ill)  an  upper  bound  n  on  the  order  of  the  plant  transfer  function, 
and  that 

Civ)  the  seros  of  the  plant  transfer  function  lie  in  the  open  left  half 

plane. 

Of  these,  conditions  (1)  and  (11)  are  quite  restrictive.  As  might  be  expected, 
the  corresponding  conditions  for  the  multivariable  systems  are  considerably  more 
stringent.  The  principal  aim  of  this  paper  are:  ( 


mi  i*  %  «rs  .  £  ^  ... 

f  ij  •’  i  **  <t 


2 


1)  to  elucidate  these  conditions  and  to  discuss  the  considerations  which 
arise  in  the  design  of  globally  stable  adaptive  controllers  for  mxm  multi-input 
multi-output  (MIMO)  systems, 

and  il)  to  examine  the  nature  of  the  prior  information  needed  for  a  complete 
solution  of  the  adaptive  control  problem  for  2x2  systems. 

In  section  IVa  it  is  shown  that  the  knowledge  of  the  relative  degree  of 
an  S1S0  system  generalizes  to  the  knowledge  of  the  Hermite  normal  form  of  the 
plant  transfer  matrix.  The  Hermite  form  plays  a  central  role  in  the  choice 
of  a  reference  model  in  the  multivariable  case.  In  section  Vb  the  definiteness 
of  a  gain  matrix  associated  with  the  plant  transfer  matrix  is  shown  to  correspond 
to  the  multivariable  version  of  condition  (il)  and  is  needed  to  generate  stable 
adaptive  laws.  Hence,  the  feasibility  of  using  adaptive  control  for  the  MIMO 
case  in  practical  situations  hinges  strongly  on  the  availability  of  prior 
information  needed  to  satisfy  the  above  conditions. 

Early  attempts  to  extend  SISO  results  to  the  MIMO  case  were  made  by 
Monopoll  and  Hsing  [4J  for  continuous  time  systems  and  by  Borlson  [5]  and  Goodwin 
et  al.  [6]  for  discrete  time  systems.  All  of  them  tacitly  assumed  that  the 
plant  transfer  matrix  can  be  diagonalized.  More  recently.  Monopoll  and  Subbarao 
[10]  have  considered  a  special  class  of  such  2x2  systems  for  practical  applica¬ 
tions.  In  [7]  Elliott  and  Wolovlch  Introduced  the  concept  of  the  Interactor 
[15]  and  later  Goodwin  and  Long  used  this  concept  to  generalize  the  results  In 
[8]  .independently  of  this  work,  recently,  Morse  [9]  discussed  the  Importance  of 
the  Hermite  normal  form  of  a  transfer  matrix  defined  over  the  principal  ideal 
domain  (PID)  of  proper  rational  functions  [11]  in  the  context  of  general  MIMO 
adaptive  control.  In  [9]  it  is  stated  that  the  Hermite  form  and  the  Interactor 
contain  equivalent  information. 
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As  mentioned  earlier,  the  Hermlte  form  (or  equivalently  the  Interactor) 

of  the  plant  transfer  matrix  and  the  sign  definiteness  of  the  high  frequency 

gain  matrix  have  to  be  known  aprlori  before  a  stable  adaptive  controller  can 

be  designed.  In  general,  as  discussed  in  section  IVa,  the  Hermlte  form  has  a 

triangular  structure.  In  the  approach  used  here,  it  is  assumed  that  the  relative 

degree  of  each  element  of  the  plant  transfer  matrix  is  known.  This  information 

Is  adequate  to  determine  whether  the  Hermlte  form  is  diagonal  or  triangular, 

and  whether.  In  the  latter  case,  it  can  be  made  generlcally  diagonal  using  a 

known  prefilter.  When  the  Hermlte  form  is  triangular,  its  off-diagonal  elements, 

in  general,  depend  upon  the  unknown  plant  parameters  and  hence  can  not  be 

specified  aprlori.  Therefore,  a  sufficient  condition  for  adaptive  control  to 

be  practically  feasible  is  that  the  Hermlte  form  be  diagonal.  Even  when  this 

condition  is  satisfied,  the  high  frequency  gain  matrix  K  can  be  either  diagonal 

P 

or  non-diagonal.  In  the  former  case  only  the  sign  of  each  diagonal  element  needs 
to  be  known  for  generating  stable  adaptive  laws.  When  is  not  diagonal,  the 

additional  prior  Information  regarding  the  definiteness  of  its  symmetric  part 
must  be  available  for  the  adaptive  control  of  the  multivariable  plant.  These 
different  cases  are  illustrated  by  considering  2x2  systems  in  detail.  It  is 
shown  there  that  all  stably  invertible  2x2  plants  can  be  generlcally  adaptively 
controlled  subject  to  the  definiteness  of  the  gain  matrix. 

Section  II  states  the  problem  of  multivariable  adaptive  control  in  a  general 
setting.  Section  III  contains  four  Important  lemmas  from  the  adaptive  control 
and  multivariable  literature.  These  are  used  extensively  in  the  subsequent  two 
sections  in  setting  up  the  reference  model  and  realizing  the  controller  in  the 
feedback  form  and  in  proving  the  stability  of  the  overall  system  for  different 
types  of  2x2  plants  in  section  V. 
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II.  Sta  tenent  of  the  Problem 

An  mxm  multi-input  multi-output  linear  time  Invariant  plant  P  is  completely 
represented  by  the  m-lnput  m-output  vector  pairs  (u(») »yp(*) J •  It  is  assumed 
that  P  can  be  modeled  by  a  rational  transfer  matrix 

W  (s)  -  Z  Cs)R"1(s)  (1) 

P  P  P 

with  {Z  (s)  and  R  (s)}  right  coprime  polynomial  matrices,  both  of  dimension 
P  P  ' 

mxm  and  Rp(s)  column-proper  (i.e.,  the  constant  matrix  [Rp]h formed  by  the  co¬ 
efficients  of  the  highest  powers  in  each  column  of  R^(s)  is  nonsingular) .  Further, 

W  (8)  is  of  full  rank  and  is  strictly  proper.  The  zeros  of  the  plant  transfer 
P 

matrix,  given  by  the  roots  of  the  polynomial  det  [Z^(s) ] ,  lie  In  the  open  left 
half  plane,  while  the  plant  poles  may  be  unstable.  The  parameters  of  Wp(s)  are 
assumed  to  be  unknown. 

A  reference  model  represents  the  behavior  expected  from  the  plant  when  the 
latter  is  augmented  with  a  suitable  differentiator  free  controller  (a  cascade 
controller  in  combination  with  linear  state  feedback) .  The  model  is  linear  time- 
invariant  and  has  a  piecewise  continuous  and  uniformly  bounded  reference  input 
vector  r(* )  and  output  vector  y^C* ) •  The  transfer  matrix,  denoted  by  W^(s) , 
is  strictly  proper  and  stable. 

The  error  between  the  plant  and  the  model  outputs  is  defined  as 

e.  Ct)  ■  y  ( t)  -  y  (t) .  C2) 

1  p  m 

The  adaptive  control  problem  is  to  determine  a  suitable  control  vector  u(*) 
such  that 

lim  ||e.(t)H  -  lia  |y  Ct)  -  y  (t)|  -  0  (3) 

t-*»  t-"°  v 


As  in  the  scalar  case,  the  solution  to  the  above  problem  can  be  divided 
into  two  parts  -  an  algebraic  part  and  an  analytic  part.  The  algebraic  part 
is  concerned  with  the  equivalence  class  of  reference  models  which  can  be  used 
as  well  as  the  structure  of  the  adaptive  controller  whose  parameters  are  to 
be  adjusted.  The  existence  of  a  solution  then  corresponds  to  the  existence 
of  a  constant  controller  parameter  matrix  such  that  condition  (3)  is  satisfied 
for  any  arbitrary  input  r(*).  The  structure  of  the  controller  also  determines 
the  uniqueness  (or  nonuniqueness)  of  the  solution. 

Once  the  existence  of  a  solution  is  established,  the  analytic  part  deals 
with  adaptive  schemes  for  updating  the  unknown  control  parameter  matrix,  so 
that  error  e^(t)  evolves  asymptotically  to  zero. 

HI  •  Mathematical  Preliminaries 

The  algebraic  and  analytic  aspects  of  the  adaptive  control  problem  are 
discussed  in  the  following  two  sections  and  use  well  known  results  in  linear 
multivariable  theory  and  stability  theory  of  dynamical  systems.  These  results 
are  presented  here  as  four  principal  lemmas  and  their  relevance  to  the  multi- 
variable  adaptive  control  problem  is  briefly  discussed.  Lemma  1  (Bezout  Identity 
for  polynomial  matrices)  is  proved  here  following  the  proof  given  in  [1]  for 
scalar  polynomials.  Proofs  of  lemmas  2  and  3  can  be  found  in  standard  texts 
on  multivariable  systems  [13,14].  A  brief  outline  of  the  proof  of  lemma  4,  which 
is  the  multivariable  version  of  the  Lln-Narendra  error  model  [17],  is  presented 
here  for  easy  reference, 
a)  Bezout  Identity 

Lemma  1:  Let  Q(s)  and  T(s)  be  mxm  right  coprime  polynomial  matrices  with 

each  column  degree  of  T(s)  strictly  less  than  the  corresponding  column  degree 

•■1 

dj  of  Q(s),  with  Q(s)  column  proper  (i.e.,  T(s)Q  (s)  is  a  strictly  proper 
transfer  matrix).  Then  mxm  polynomial  matrices  P(s)  and  R(s),  each  having 


highest  degree  (v-1),  exist  such  that  P(s)Q(s)  +  R(s)T(s)  can  be  made  equal  to 
any  arbitrary  mxm  polynomial  matrix  M(s)  with  each  column  degree  less  than  or 
equal  to  (dj  +  v-1) >  where  v  is  the  observability  index  of  the  minimal  transfer 
matrix  T(s)Q  1(s) . 

Proof:  Since  Q(s)  and  T(s)  are  right  coprime  polynomial  matrices,  there  exist 
polynomial  matrices  A(s)  and  B(s)  [12]  such  that 

A(s)Q(s)  +  B(s)T(s)  -  I. 

Let  M(s)  be  an  mxm  arbitrary  polynomial  matrix  with  column  degree  $  (d^  +  v-1) . 
Then 

M(8)A(s)Q(8)  +  M(s)B(s)T(s)  -M(s).  (4! 

The  right  coprlme  factorization  T(s)q“1(«)  can  also  be  expressed  by  a  left 
coprime  factorization  E**1(s)F(s)  with  E(s)  row  proper  and  each  row  degree  of 
F(s)  strictly  less  than  the  corresponding  row  degree  of  E(s) .  The  highest 
degree  of  ECs)  is  v,  the  observability  index  of  the  minimal  transfer  matrix 
[14].  Hence, 

T(s)Q‘1(s)  -  E_1(s)F(8)  (5 

(4)  and  (5)  can  be  represented  as  a  composite  matrix  equation 


M(s)A(s) 

M(s)B(s5|  [~q 

(8)1  [m(s‘) 

■1 

_  F(s) 

-E(s)  J  Lt 

Cs]J  Lo. 

By  elementary  column  operations,  the  ahove  matrix  equation  can  be  reduced  to 


such  that  every  column  degree  of  R(s)  Is  strictly  less  than  the  corresponding 
column  degree  of  E(s) .  Since  the  highest  degree  of  E(s)  is  v,  the  highest  degree 
of  R(s)  can  be  at  most  (v-1) .  Further,  since  Q(s)  is  column  proper  and  each 
column  of  M(s)  has  degree  less  than  or  equal  to  (d^  +  v  -  1) ,  the  polynomial 
matrix  P(s)  can  have  degree  at  most  (v-1) . 

This  lemma  is  used  to  establish  the  existence  of  a  controller  structure  so 
that  the  transfer  matrix  of  the  plant  together  with  the  controller  is  Identical 
to  that  of  the  model, 
b)  Decoupling  by  State  Feedback 

Let  G(8)  be  a  nonsingular  mxm  strictly  proper  transfer  matrix.  Let  d^ 
denote  the  minimum  relative  degree  in  the  ith  rov  of  G(s) ,  i.e.,  <i,  •  min  (degree 
difference  in  s  of  the  denominator  and  the  numerator  of  each  entry  of  the  1th 
rov  of  G(8>)  -  1.  Let  (lxm)  constant  row  vector  E^  be  defined  as 

A  di+l 

E  -  lim  s  G(s)  . 

1  s-*- 

It  is  known  {13J  that  d^  and  E^  are  invariant  under  linear  state  feedback. 

Lenma_2:  Let  G(s)  and  E^,  1  “  l,...,m  be  defined  as  above.  G(s)  can  be  decoupled 
by  linear  state  feedback  if  and  only  if  the  constant  matrix 


is  nonsingular. 

The  entries  in  the  matrix  E  are  the  high  frequency  scalar  gains  associated 
with  scalar  transfer  functions  of  minimum  relative  degree  in  each  row  of  the 
transfer  matrix.  This  lemma  is  used  in  section  Va  to  specify  the  model  from  a 
knowledge  of  the  relative  degree  of  each  entry  in  the  plant  transfer  matrix. 
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c)  Number  of  Inputs  and  Outputs 

Lemma  3:  A  transfer  matrix  G(s)  is  output  function  controllable  if  and  only  if 
it  has  linearly  independent  rows  over  the  field  of  rational  functions,  i.«.,  the 
rank  of  G(s)  is  equal  to  the  number  of  outputs. 

In  model  reference  adaptive  control,  for  the  plant  output  to  follow  the 
model  output  asymptotically,  the  plant  transfer  matrix  must  be  output  function 
controllable.  If  such  a  transfer  matrix  has  more  columns  than  rows,  then  the 
number  of  inputs  in  excess  of  the  number  of  outputs  of  the  plant  can  be  set  to 
an  arbitrary  value,  and  in  particular,  to  zero.  This  pertains  to  the  columns 
of  zeros  in  the  Hermite  normal  form  [cf.  11]  of  a  rectangular  transfer  matrix. 
Alternatively,  this  is  equivalent  to  selecting  inputs  corresponding  to  linearly 
independent  columns  of  the  transfer  matrix,  which  in  turn  is  reflected  in  the 
nonzero  columns  of  its  Hermite  form.  Hence,  in  general,  it  is  sufficient  to 
consider  square  transfer  matrices  with  the  same  number  of  Inputs  and  outputs. 

d)  Multivariable  Error  Model  Prototype  3 
Lemma  4 

Given  a  stable  m-lnput,  m-output  n-dimenslonal  minimal  triple  (C,A,B) ,  two 
symmetric  positive  definite  matrices  T  and  and  u»(t)  :  [0,»>  ■+  RP  whose  elements 
are  piecewise  continuous  functions,  the  equilibrium  state  of  the  set  of  differ¬ 
ential  equations 

e(t)  -  Ae(t)  +  Bv(t) 
ex<t)  -  Ce(t) 

v(t)  -  *(t)a>(t)  -  WT(t)rim(t)e1Ct)  (7) 

♦(t)  -  -rei<t)u>T(t)  (8) 

is  stable  if  the  transfer  matrix  T(s)  -  C(sl-A)  *B  is  strictly  positive  real  (SPR) . 


Proof : 


The  proof  follows  directly  by  using  a  Lyapunov  function  candidate 

V(e,*)  ^  eT(t)Pe(t)  +  trace  [#T(t)r  1*Ct)  ]  with  P  -  PT  >  0  (9) 

and  the  matrix  version  of  the  Kalman-Yacubovich  Lemma  [18].  This  yields 

V(e,*)  -  -eT(t)[GGT  +  eN]e(t)  -  2e  T(t)e  (t)»T(t)r  »(t)  <  0  (10) 

T 

for  some  G  and  H  *  M  >0  matrices  and  e  >  0. 

It  follows  that  the  error  model  is  uniformly  stable  and  e(t)  (hence  e^(t)) 
and  ♦( t)  are  uniformly  bounded  for  all  finite  initial  conditions.  Further,  from 
(g)  and  (ia)it  follows  that  e(*)  and  ♦(•)et  •  If  in  addition  it  is  assumed  that 
u)(t)  and  u(t)  are  uniformly  bounded,  it  can  be  concluded  that  lim  e(t)  “  0 
and  lim  *(t)  ■  0.  However,  very  little  can  be  said  about  the  convergence  of 

*(•)  to  a  constant  matrix. 

Remark:  If  T(s)  ■  1  in  the  above  lemma,  the  third  error  model  (eqns.  7,8)  degener¬ 
ates  into  the  first  error  model.  The  equations  describing  such  a  model  and  the 
corresponding  adaptive  equations  may  be  expressed  as  follows 


*(t)w(t) 

-  eL(t) 

(11) 

e 

re  (t)<*)T(t) 

T 

r  -  r  >  o 

*(t)  -  - 

1 

- T - 

T 

1-ho  (t)T^u(t) 

F1  “  ri  >  °‘ 

(12) 

e 

Using  similar  arguments,  it  can  be  shown  that  when  w(t)  and  ui(t)  are  uniformly 
bounded,  lim  e.(t)  -  0  and  lim  *(t)  ■  0. 


IV.  Structure  of  the  Adaptive  S ye ten:  General  Case 
a)  Hermite  Normal  Form 

In  a  model  reference  adaptive  system,  the  output  of  the  plant  is  required  to 

follow  the  output  of  a  reference  model.  The  plant  together  with  the  controller, 

whose  parameters  are  adjusted,  asymptotically  approaches  a  linear  time-invariant 

system.  The  transfer  matrix  of  the  latter  should  therefore  be  identical  to  that 

of  the  reference  model  for  perfect  model  following,  if  the  class  of  inputs  is 

sufficiently  general.  An  important  question  that  has  to  be  resolved  in  the  initial 

stages  of  design  is  the  choice  of  the  model  transfer  matrix,  i.e.,  the  class  of 

rational  transfer  matrices  from  which  the  model  transfer  matrix  should  be  selected. 

It  is  in  this  context  that  the  Hermite  normal  form  is  found  to  be  Important. 

The  principal  idea  here  is  that  the  set  of  all  stable  reference  models  can  be 

generated  by  the  stable  Hermite  form  of  the  plant  transfer  matrix. 

The  following  concepts  from  linear  systems  theory  are  found  to  be  relevant 

for  a  discussion  of  the  Hermite  form. 

The  set  of  all  proper  rational  transfer  function  Rp(s)  is  a  principal 

ideal  domain  (PID)  .  A  matrix  with  elements  in  R  (s)  is  invertible  if  and  only 

P 

if  its  determinant  is  a  unit  in  R  (s) .  A  unimodular  matrix  is  an  Invertible 

P 

proper  transfer  matrix  whose  inverse  is  also  a  proper  transfer  matrix.  If  the 

'relative  degree'  (degree  of  the  denominator  minus  the  degree  of  the  numerator) 

of  each  rational  transfer  function  is  taken  as  the  'degree'  of  each  element,  a 

division  rule  can  be  established,  making  this  a  Euclidean  domain  [16]. 

Two  transfer  matrices  T^s)  and  T2(s)  over  Rp(s)  are  said  to  be  dynamically 

equivalent  if  and  only  if  there  exists  a  unimodular  matrix  C(s)  over  R  (s)  such 

P 

that  T^Cs)  -  T.(s)C(s)  .  The  Hermite  normal  form  for  nonsingular  matrices  over 
PID  [12],  obtained  by  performing  elementary  column  operations  on  the  matrix, 
represents  the  canonical  form  within  each  equivalence  class.  This  together  with 
the  rank  of  the  transfer  matrix  represents  the  complete  set  of  invariants  within 
each  class  [11]. 
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The  Hecnite  form  of  an  mxm  matrix  T(s)  over  IR^Cs)  la  a  lower  triangular 
rational  matrix  of  the  form 


H(s)  - 


where  h, ,(s)  -  { ,  nj(  <  n  },  is  proper  and  n  and  n  are  positive 

y  ■>„  id  i  1  « 

Integers.  ir(s)  is  any  monic  polynomial  of  degree  1.  The  choice  of  ir(s)  is 
immaterial,  but  once  it  is  chosen  the  Hermite  form  H(s)  is  unique.  If  H^(s) 
is  n^-Hermite  form  and  H^(s)  is  ^"Hermite  form,  then  either  can  be  obtained 
from  the  other  by  elementary  column  operations.  In  other  words,  H^(s)  and  H^Cs) 
are  dynamically  equivalent.  However,  if  the  roots  of  tt(b)  lie  in  the  open  left 
half  plane,  H(s)  corresponds  to  a  stable  Hermite  form. 

From  the  foregoing  discussion  it  follows  that  every  plant  transfer  matrix 
Wp(s)  generates  a  class  C  of  dynamically  equivalent  models.  The  set  of  all 
stable  reference  models  which  the  plant  transfer  matrix  can  follow,  denoted  by 
(l/,  is  a  proper  subset  of  C.  Hence,  by  the  preceding  considerations,  a  stable 
Hermite  form  of  the  plant  transfer  matrix  itself  can  be  chosen  as  a  reference 
model.  The  entire  class  W  of  the  stable  reference  models  can  then  be  generated 
by  pos multiply ing  the  Hermite  form  by  a  known  and  fixed  dynamic  controller. 

The  Importance  of  the  Hermite  form  H(s)  e  W  lies  in  the  fact  that  it  can  be 
determined  a  priori  directly  from  reasonable  Information  about  the  plant  trans¬ 
fer  matrix  V  Cs) .  In  algebraic  terms,  the  Hermite  form  represents  a  basis  for 
P 

the  free  module  over  R  (a)  spanned  by  the  columns  of  the  plant  transfer  matrix. 

P 

Dynamic  equivalence  then  Implies  that  the  free  modules  spanned  by  the  columns 


of  the  plant  and  model  transfer  matrices  are  the  same. 


For  SISO  systems  the  Hermite  form  is  simply  — *  >  where  ir  is  as  defined 

before  and  n  is  the  relative  degree  of  the  transfer  function.  The  unlmodular 

form  is  lust  a  unit  in  K  (s) .  The  invariance  of  the  relative  degree  and  the 

P 

realization  of  the  controller  in  the  feedback  and  the  feedforward  form  are  well 

known  in  adaptive  control  literature.  In  that  sense  deg[det  H(s)]  represents 

the  multivariable  analog  of  the  relative  degree  (deg[det  R  (s) ]  minus  deg[det  Z  (s)]) 

P  P 

It  is  shown  in  the  next  subsection  how  the  multivariable  controller  can  be  realized, 
b)  Controller  Structure 

From  the  discussion  in  section  IVa  it  follows  that  the  reference  model  trans¬ 
fer  matrix  W  (s)  is  dynamically  equivalent  to  the  plant  transfer  matrix  W  (s),  i.e., 
m  P 

W  (s)Q(s)  -  W  (s)  (14) 

p  m 

for  some  unlmodular  matrix  Q(s) .  If  H(s)  is  the  Hermite  form  in  the  equivalence 
class,  then 


wp(s)Qp(s)  -  H(s)  -  Wm(a)Q^(«) 


where  Qp(s)  and  (^(s)  are  unlmodular  matrices.  From  the  above  two  equations  it 
readily  follows  that 


w  (s)Q  (s)(T  Cs)  -  W  (s)Q(s)  -  W  (s). 
ppm  p  m 


The  basic  structure  of  the  adaptive  system  is  shown  in  Fig.  1.  The  feed¬ 
forward  controller  Q  \s)  is  fixed  and  known.  The  controller  Q  (s)  can  be 

m  p 

realized  as  shown  below. 

Let  Qp(s)  be  factorized  as  Rp(s)PF1(s)K0  ■  Rp(s)  (K^P^s))  1  such  that 
3  JR  (s)l  -  3JP„(s)] 


and 


(18) 


11m  Q  (s)  -  11m  R  (s)Pp*(s)K  -  K  . 

8^=0  P  g^o,  P  r  0  0 

3^^ [• ]  denotes  the  column  degree  of  a  polynomial  matrix  and  [ • ]^  represents  the 

constant  matrix  formed  by  taking  the  coefficients  of  the  highest  powers  in  each 

column  of  a  matrix.  From  eqn.  (15) 

tH"1(g)w  (s)]Q  (s)  -  I  .  (19) 

P  P 

Since  Qp(s)  is  unimodular,  [H  ^(s)Wp(s)]  is  also  unimodular.  This  implies 

lim  [h”1(s)W  (s)  ]  -  lim  [(H_1(s)Z  (s))  (r”1^))  ]  -  K  (20) 

g-H*>  P  8~*CO  P  P  P 

where  Rp  is  a  constant  nonsingular  matrix,  and 

3  .[H"1(s)Z  (a)]  -  »  AR  (•)]. 
cj  p  cj  p 

is  the  high  frequency  gain  matrix  of  the  transfer  matrix.  In  view  of  (17) 
and  (19) ,  let 

K"lpFCa)  -H~1(g)Zp(s)  (21) 

such  that 


3cj[H'  Cs)Zp(.)]  -  3cj[Rp(»)]  -  3cJ[Pf(8)]  -  dj. 


From  (19) 


11®  (H_1(s)Z  (s)R"1(s))Q  (s) 
«-»•  P  P  P 


-  K  *K.  -  I 
P  0 


lim  (H  *(s)Z  (s))R  *(s)*lim  Q  (s) 
S-*»  P  P  g-H»  ? 


Hence,  K.  -  K  . 

0  p 

Then  (21)  becomes 


P_(s)  -  K"1H'1(s)Z  (s)  . 

F  p  p 


The  overall  transfer  matrix  Is  given  by 


[Z  Cs)R'1(s)]lR  (sjP^CsJKJ  -  Z  (s) [Z_1H(s)K  ]Kn  -  H(s)  (22) 

pp  p  r  U  p  p  p  U 

The  controller  Q  (s)  -  R  (slP^CsjK.  consists  of  gain  matrix  K_  In  the  forward 
p  p  F  U  0 

path  and  two  auxiliary  signal  generators  F^  and  In  the  feedback  path.  F 
contains  a  system  with  transfer  matrix  W^(8)  -  N  *(s)C(s)  and  F^  contains  a 
system  with  transfer  matrix  W^s)  -  N  *(a)D(s)  +  BQ.  For  constant  values  of 
parameters,  the  overall  transfer  matrix  from  v(t)  to  y^(t)  is  W(s) ,  where 

W(s)  -  Zp(s)R"1(s)Rp(s)[(N(s>  +  C(s))Rp(s)  +  (D(s)  +  D0N(s))Zp(8)  f1 

•N(s) *Kq  .  (23) 

By  Lemma  1,  given  a  polynomial  matrix  P  (s)  of  column  degree  d  and  an  arbitrary 

r  J 

polynomial  matrix  N(s)  with  3rjlN(s)]  -  v-1  such  that  3^ [N(s) *Pp(s) ]  “  d j  +  v"1» 
polynomial  matrices  C(s),D(s)  and  a  constant  matrix  Dq  of  appropriate  degrees  can 
be  determined  from  the  following  Bezout  identity. 

!N(s)  +  C(s) ]R  (8)  +  [D(s)  +  DnN(s) ]Z  (s)  -  N(s)P_(s) . 

p  up  t 

v  is  the  observability  index  of  the  plant  transfer  matrix  which  is  assumed  to  be 
known.* 

The  arbitrary  mxm  polynomial  matrix  N(s)  is  chosen  such  that  it  is  row 
proper  and  it  commutes  with  both  C(s)  and  D(s) .  Hence,  N(s)  Is  chosen  as  diagonal 
[n(s)J  where  n(s)  is  an  arbitrary  monic  Hurwitz  polynomial  of  degree  (v-1). 

F^  generates  a  set  of  auxiliary  signals 

•  Jet)  ■  2«T(t),...,s°uT(t)]  i  -  2,...,v  (24) 

and  contains  (v-1)  mxm  parameter  matrices  C^d  -  l,...,v-l).  Y^  generates  a  set 

*  An  upper  bound  on  v  is  sufficient,  but  then  the  above  equation  may  not  have 
a  unique  solution. 


of  auxiliary  signals 


“^(t)  "  nboIsV”lypT(t} . 8°ypCt)]  i-v+l,...,2v  (25) 

and  contains  v  mm  parameter  matrices  Dj(i  -  0,1,..  .  Together  with 

^(t)  -  y(t),  output  of  the  feed  forward  controller,  and  Kq  these  constitute 
2v  m-vector  signals  and  2v(mxm)  matrices  of  adjustable  parameters  of  the  con¬ 
troller  denoted  by  the  elements  of  a  parameter  matrix 

9(t)  -  tK0(t);C1(t),...,Cv_1(t);D0(t),D1(t),...,Dv_1(t)]  (26) 

T 

The  control  Input  u(t)  to  the  plant  is  given  by  0(t)u)( t)  -  e^(t)u>  (t) r^oi(t) 

T  * 

where  >0.  For  6(t)  ■  0  ,  the  desired  feedback  controller  parameter 

matrix,  the  plant  together  with  the  feedback  controller  yields  its  Hermite  form. 
The  parameter  error,  matrix  +(t)  is  defined  as 

♦(t)  -  0(t)  -  0* 

and  this  gives  rise  to  the  error  transfer  matrix, 
c)  Error  Equation 

The  error  transfer  matrix  representing  the  error  equation  between  the  plant 
and  its  Hermite  form  is  given  by 

W  (a)  -  H(s)K  (27) 

•  P 

Hence,  the  error  model  can  be  represented  as  shown  in  Fig.  2  with 
T 

♦( t)u(t)  -  e1(t)o>  (t)Tju(t)  as  the  input  to  a  system  with  transfer  matrix  We(s) . 
It  is  worth  noting  that  the  signals  u^(t)(l  -  l,...,2v)  are  derived  from  v(t), 
output  vector  of  a  known  and  fixed  feedforward  compensator  Qj|^(s)  ,u(t) ,  the 
plant  control  Input  vector  and  y  (t),  the  output  vector  of  the  plant. 


H(s) 


The  Error  Model 
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V.  Adaptive  Control  Problem:  2x2 

a.  In  this  subsection  the  concepts  discussed  in  section  IV  are  applied  to  the 
adaptive  control  of  2x2  systems.  It  Is  assumed  that  the  relative  degree  of  each 
of  the  four  scalar  transfer  functions  In  the  plant  transfer  matrix  is  known  and 
that  adequate  information  regarding  the  matrix  E  (cf .  eqn.  (6))  is  available  so 
that  the  nature  of  the  Hermite  form  can  be  concluded  using  lemma  2.  The  differ¬ 
ent  situations  that  can  arise  are  discussed  and  it  is  shown  through  examples 
that  adaptive  control  is,  in  general,  feasible  only  when  the  Hermite  form  is 
diagonal.  In  other  words,  the  plant  transfer  matrix  is  decoupled  using  the 
controller  structure  described  in  the  previous  section.  In  such  cases,  the 
exponents  n^  along  the  diagonal  in  H(s)  represent  the  minimum  relative  degree  in 
the  corresponding  rows  of  the  plant  transfer  matrix.  In  the  discrete  case  a 
diagonal  transfer  matrix  implies  that  each  output  is  affected  by  an  independent 
input  Cor  inputs)  with  minimum  delay.  The  nonzero  entries  in  the  high  frequency 
gain  matrix  in  (eqn.  20)  are  the  scalar  gains  associated  with  the  transfer 
functions  of  minimum  relative  degree  in  each  row  (cf.  lemma  2). 

Throughout  this  section  the  term  'controller'  will  refer  to  the  feedback 
controller  whose  parameters  are  to  be  adjusted  and  the  term  'model'  will  be 
used  interchangeably  with  the  Hermite  form  of  the  plant  transfer  matrix. 


represent  the  plant  transfer  matrix  in  which  and  (i  -  1,2, 3,4)  are  unknown 
parameters.  It  is  further  known  that  the  zeros  of  the  plant  transfer  matrix  lie 
in  the  open  left  half  of  the  complex  plane.  In  the  following  three  cases,  the 
prior  information  that  is  assumed  regarding  the  unknown  parameters  is  successively 


greater. 
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Case  (1)  Let 

In  this  case  the  model  (Hermite  form)  HCs)  is  diagonal  and  the  high  frequency 

gain  matrix  K  Is  the  same  as  E  in  eqn.  6  and 
P 


The  singularity  of  the  matrix  E  for  W^(s)  implies  that  the  model  will  be 
triangular.  The  Hermite  form  H(s)  and  the  corresponding  high  frequency  gain 
matrix  are  respectively: 


It  is  clear  that  considerable  Information  regarding  the  plant  parameters  will 
be  needed  to  set  up  the  model  as  well  as  to  generate  the  adaptive  laws. 

Case  (ill)  Even  greater  knowledge  of  the  plant  parameters  will  be  needed  for 
adaptive  control  when  k^k^  ■  k^k^  and  ■  Oj  +  resulting  in 


1 

s+a 

0 

“l 

k2 

H(s)  - 

•j^ls+a+Cdj-aj)  ] 

1 

a  >  0;  K  • 
P 

k3<“l-V 

k4(o2"a4) ' 

(s+a) 1 

(e+a)3 

(a-o3> 

(s-04) 

Hence,  even  in  the  simple  case  described  in  this  example,  adaptive  control  may  be 
practically  feasible  (in  terms  of  the  prior  information  needed)  only  in  case  (1) 
where  H(s)  is  diagonal  and  K  has  a  simple  structure. 


■  V 


Example  2  In  this  cam  the  transfer  matrix 


y*> 


(»4-04) 


and  has  transfer  functions  of  minimal  relative  degree  in  both  rows  in  the  first 
coltmm.  Once  again  the  matrix  E  is  singular  so  that  the  model  H(s)  is  triangular 
and  the  high  frequency  gain  matrix  Kp  is  no  longer  the  same  as  E. 


H(s)  - 


3 

(s+a) 


(s+a) 


a  >  0  ;  K  - 
P 


V1 


For  the  same  reasons  as  before,  considerable  prior  information  regarding  the  unknown 
plant  parameters  will  be  needed  in  this  case  also  for  the  implementation  of  an 
adaptive  controller. 

The  above  two  simple  examples  illustrate  the  two  different  ways  in  which  the 
Hermlte  form  of  a  plant  transfer  matrix  may  become  triangular.  By  completely 
classifying  (2x2)  transfer  matrices  in  terms  of  their  Hermlte  forms  we  can  es tab- 
llsh  the  prior  Information  needed  to  control  tbem  adaptively. 

Let  the  transfer  matrix  of  a  2x2  plant  be  represented  by 


W  (s)  - 
P 


“l  °2 


k  !* 


where  kA  (i  -  1,2, 3,4)  are  the  high  frequency  scalar  gains  and  (i  •  1,2, 3,4) 
are  the  relative  degrees  of  the  four  transfer  functions.  Using  lemma  2,  and 


20 


knowledge  of  ,  the  models  and  the  gain  matrices  corresponding  to  four  differ¬ 
ent  classes  of  all  2x2  plants  are  delineated  in  Table  1.  The  number  of  elements 

*  *  *  * 

in  each  class  is  shown  in  column  1  (for  example  n^  <  n2  ’  °3  >  n4  and 

*  *  *  * 

n,  >  n~  ,  n  <  n.  are  the  two  elements  of  class  I) .  The  structure  of  only 
1*34 

one  typical  element  in  each  class  is  given  in  column  2  and  every  element  of 

the  class  can  be  reduced  to  the  same  form  by  relabling  the  rows  and/or  columns 

of  W  (s)  . 

P 

In  class  1  both  the  model  H(s)  and  the  gain  matrix  K  are  diagonal.  Adaptive 

P 

control  is  possible  in  this  case  if  the  sign  of  the  elements  of  are  known.  In 
class  II,  which  contains  four  elements,  H(s)  is  diagonal  but  Kp  is  triangular. 


Condition  on  the  Structure  of 


rel.  deg. 

*  * 

I  <"2 

2  elements  *  * 

“3  >n4 


*  * 
“l  <n2 


4  elements 


n3  ’  °4 


*  * 

III  n.  “  n„ 

1  2 


1  element  * 


n2  °4 


*  * 

IV  “l  <n2 
2  elements  *  * 

n3  <n4 


Model,  H(s) 


(s+a) 


(s+a) 


(s+a) 


(»+a) 


n_  •n. 
3  4 


(s+a) 


n  -n_ 
1  2 


(s+a) 
Non  diagonal 


n„  +n/ 
3  4 


Structure  of 
gain  matrix  K 


\\  *  k2k 


Not  obvious 


Table  I 


A  sufficient  condition  for  adaptive  control  is  that  K  be  sign  definite  (shown 

P 

in  the  next  subsection)  and  hence  considerably  more  prior  information  is  needed 

in  this  case.  In  class  111,  which  has  one  element,  the  minimum  relative  degree 

in  each  row  occurs  in  both  columns.  If  k, k  i  k„k  and  the  matrix  K  is  sign 

1  4  2  3  P 

definite,  the  plant  can  be  adaptively  controlled. 

In  class  IV  the  transfer  functions  with  minimum  relative  degree  (in  each 
row)  occurs  in  the  same  column  so  that  H(s)  is  triangular  in  structure.  The 
corresponding  gain  matrix  also  depends  in  general  u.on  the  unknown  plant 
parameters  in  a  complex  way.  Hence,  for  this  class  of  transfer  matrices, 
adequate  prior  information  is  generally  not  available  to  enable  adaptive  con¬ 
trollers  to  be  designed  directly.  However,  as  shown  in  example  3,  using 
dynamic  compensation,  the  Hermite  form  of  the  modified  transfer  matrix  may  be 
made  diagonal.  The  following  simple  example  illustrates  how  this  can  be  achieved 
Example  3:  In  example  2  it  was  shown  that  when 


the  resulting  H(s)  would  be  triangular  and  contain  the  unknown  plant  parameters. 

-L_  0 

If  such  a  plant  is  augmented  by  a  matrix  D(s)  ■  *+8 


,  B  >  0,  then  the 


L_° 


modified  plant  transfer  matrix 


W  vs)  -  W  (s) D(s)  - 
p  P 


Tif;3^»+6)  (rta.)2 


has  a  diagonal  Hermite  form  n(s)  where 


A  sufficient  condition  for  adaptive  control  is  that  K  be  sign  definite  (shown 

P 

in  the  next  subsection)  and  hence  considerably  more  prior  information  is  needed 

in  this  case.  In  class  III,  which  has  one  element,  the  minimum  relative  degree 

in  each  row  occurs  In  both  columns.  If  k.k  +  k  k  and  the  matrix  K  is  sign 

14  <3  P 

definite,  the  plant  can  be  adaptively  controlled. 

In  class  IV  the  transfer  functions  with  minimum  relative  degree  (in  each 
row)  occurs  in  the  same  column  so  that  H(s)  is  triangular  in  structure.  The 
corresponding  gain  matrix  also  depends  in  general  upon  the  unknown  plant 
parameters  in  a  complex  way.  Hence,  for  this  class  of  transfer  matrices, 
adequate  prior  information  is  generally  not  available  to  enable  adaptive  con¬ 
trollers  to  be  designed  directly.  However,  as  shown  in  example  3,  using 
dynamic  compensation,  the  Hermite  form  of  the  modified  transfer  matrix  may  be 
made  diagonal.  The  following  simple  example  illustrates  how  this  can  be  achieved. 
Example  3:  In  example  2  it  was  shown  that  when 


the  resulting  H(s)  would  be  triangular  and  contain  the  unknown  plant  parameters. 

. _ 0 

If  such  a  plant  is  augmented  by  a  matrix  D(s) 


*  x_ 

s+ 


s+8 

u> 


modified  plant  transfer  matrix 


W  vs)  -  W  (s)D(e)  - 
P  P 


(s+OjHs+B) 


has  a  diagonal  Hermite  form  H  (s)  where 


,  3  >  0,  then  the 


(3+02)' 


TWSPT,+6)  (rK.,)2 


A  sufficient  condition  for  adaptive  control  is  that  Kp  be  sign  definite  (shown 
in  the  next  subsection)  and  hence  considerably  more  prior  information  is  needed 
in  this  case.  In  class  III,  which  has  one  element,  the  minimum  relative  degree 
in  each  row  occurs  in  both  columns.  If  and  ths  matrix  Kp  is  sign 

definite,  the  plant  can  be  adaptively  controlled. 

In  class  IV  the  transfer  functions  with  minimum  relative  degree  (in  each 
row)  occurs  in  the  same  column  so  that  H(s)  is  triangular  in  structure.  The 
corresponding  gain  matrix  Kp  also  depends  in  general  upon  the  unknown  plant 
parameters  in  a  complex  way.  Hence,  for  this  class  of  transfer  matrices, 
adequate  prior  Information  is  generally  not  available  to  enable  adaptive  con¬ 
trollers  to  be  designed  directly.  However,  as  shown  in  axample  3,  using 
dynamic  compensation,  the  Hermlte  form  of  the  modified  transfer  matrix  may  be 
made  diagonal.  The  following  simple  example  illustrates  how  this  can  be  achieved 


,  P  >  0,  then  the 


the  resulting  H(s)  would  be  triangular  and  contain  the  unsown  plant  parameters. 

pL_ 

If  such  a  plant  is  augmented  by  a  matrix  D(s)  -  s+P 

L<> 

modified  plant  transfer  matrix 

~~  1 


W  A(s)  -  W  (s) D(s) 
P  P 


TS535T,+6)  (rh,,)2 


has  a  diagonal  Hermlte  form  HA(s)  where 


H1(8> 


(s+a) ‘ 


(s+a)_ 

The  corresponding  high  frequency  gain  matrix  K  is 

P 
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a  >  0 


0 


which  is 


the  same  as  the  matrix  E  for  W  (s). 

P 

Using  the  above  approach  we  now  show  that  all  transfer  matrices  W  (s)  which 

P 

correspond  to  class  IV  can  be  reduced  to  one  of  the  previous  classes  using  dynamic 
prefilters.  The  three  distinct  cases  which  can  arise  are  treated  separately 

below  and  the  minimal  order  diagonal  compensators  are  specified  in  each  case. 

*  *  *  * 

1)  nx  +  n4  <  n2  +  n3 

Using  a  dynamic  precompensator  with  transfer  matrix  D(s)  where 


D(s) 


Cn,  *-n  *) 


u  >  0 


the  Hermite  form  H(s)  and  the  gain  matrix  K  of  the  modified  system  can  be  obtained 

P 


as 


H^s) 


- X - X - W~ 

(nl  “n3  ) 


•  V 


,.****  . 

Cii)  ^  *  n2  +  n3  and  k^  P  k^ 


If  u^(t)  is  filtered  by 


- —H* - *  and  n,(t)  by  .1,  the  new  Hermite  form  H  (s) 

in.  -n«  )  z 

( s+ci)  *  3 


and  gain  matrix  Kp  are  given  by  the  following. 


H1(s) 


1 

— X - x 

n2  -n1 


n3  "n4| 


,  K. 


0 


(iii)  n.  +  n4  >  n2  +  n3 

Filter  u,  (t)  by  - - — * - *  and  u,Ct)  by  1  to  obtain 

J.  n«  -n.  i 

(s-hn)  1 


H1(s) 


*  *  * 
n„  +n0  -n. 
.321 


k  k 
1  2 


Remark  In  [19],  using  a  geometric  approach,  the  generic  decoupling  of  a  plant 
represented  by  the  triple  (C,A,B)  using  only  state  feedback  is  discussed.  However, 
since  the  mapping  from  the  parameter  space  of  (C,A,B)  to  the  parameter  space  of 
T(s),  the  transfer  matrix,  is  not  bicontinuous,  properties  true  in  (C,A,B)  space 
may  not  be  true  in  the  parameter  space  of  T(s) . 

Since  the  construction  of  the  Hermite  form  or  the  interactor  is  done  in  the 
parameter  space  of  the  transfer  matrix,  its  generic  decoupling  can  not  be  concluded 
on  the  above  basis.  The  problem  arises  because  the  matrix  E(Eqn.  6'  is  not  a 
continuous  function  of  its  arguments,  namely  the  minimum  relative  degr^*  in  each 
row  of  the  transfer  matrix. 

Hence,  the  point  made  in  [8]  that  the  dlagonality  of  the  interactor  matrix 
(equivalently  the  Hermite  normal  form)  is  a  generic  property  does  not  apply  here. 

In  fact,  as  shown  above,  almost  one  third  of  all  2x2  transfer  matrices  have  non¬ 
diagonal  Hermite  forms, 
b)  Adjustment  Laws: 

In  this  subsection  adaptive  laws  for  the  adjustment  of  the  controller 
parameters  are  developed  for  all  2x2  systems  categorized  in  table  1.  Since,  by 
introducing  known  prefilters,  transfer  matrices  which  belong  to  class  IV  can  be 
reduced  to  one  of  classes  I-II1,  only  the  latter  are  considered  here.  The  analysis 


is  limited  to  error  models  which  arise  in  the  various  cases.  Only  a  brief 


description  of  the  corresponding  controller  structures  is  given  since  they  are 
merely  multivariable  extensions  of  those  used  in  SISO  systems  discussed  extensively 
in  [1]  and  [2]  and  modifications  of  the  basic  structure  described  in  the  previous 
section.  The  proof  of  stability  of  the  overall  control  loop  for  SISO  systems 
in  [2]  can  also  be  directly  extended  to  these  multivariable  systems  and  hence  is  not 
discussed  here. 

The  complexity  of  the  controller  used  in  the  various  cases  depends  upon  the 

amount  of  prior  information  available  regarding  the  gain  matrix  Kp  as  well  as  the 

Hermite  form  H  (s) .  The  simplest  cases  occur  when  K  is  known  and  H  (8)  is  strictly 
P  P  P 

positive  real.  In  the  following,  the  various  cases  are  arranged  in  increasing  order 
of  generality. 

(i)  Kp  Known:  If  is  known,  considerable  simplification  is  achieved  by  including 
a  fixed  gain  matrix  in  the  control  loop.  The  modified  plant  transfer  matrix 
then  has  a  high  frequency  gain  matrix  which  is  unity  (i.e.  KQ  =  I) .  The  adaptive 
controller  contains  2x(4v-2)  parameters  which  are  the  elements  of  a  parameter 
error  matrix  $(t)  and  the  corresponding  (4v-2)  dimensional  vector  signal  is  denoted 
by  o)(t)  . 

a)  H(s)  SPR  (strictly  positive  real)*: 

When  H(s)  is  SPR  the  parameter  error  matrix  is  updated  according  to  the  law 


»(t)  -  0(t)  -  -re^(t)w(t) 


r  -  r  >o 


By  lemma  4  the  output  error  e^(t)  and  *(t)  are  bounded.  Since  the  output  of  the 
model  is  bounded,  this  ensures  (as  in  the  SISO  case)  that  the  state  of  the  plant 


and  hence  w(t)  is  bounded.  Hence  e^(t)  -»■  0  and  #(t)  +  0  as  t  +  •.  It  is 

-T 

worth  pointing  out  that  for  this  case  the  feedback  term  e^ui  (t)rco(t)  is  not  needed 
to  prove  global  stability  of  the  adaptive  loop. 


*  A  diagonal  matrix  of  rational  functions  is  strictly  positive  real  if  and  only 
if  each  diagonal  element  is  strictly  positive  real. 
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b)  H(s)  not  SPR 

If  H(s)  is  not  SPR,  auxiliary  signals  have  to  be  generated  and  added  to  the 

model  output  to  avoid  differentiation.*  The  controller  of  Fig.  1  is  modified  as 

shown  in  Fig.  3.  The  simple  error  model  of  Fig.  2.  changes  to  the  augmented 

error  model  shown  in  Fig.  4. 

T  T 

Let  8^  (t)  and  $  (t)  denote  the  ith  row  in  0(t)  and  $(t)  matrices  respectively. 

Let  hj. (s)  and  t*ie  tw0  non-zero  diagonal  elements  of  H(s)  .  In  construct¬ 

ing  the  auxiliary  signals  the  matrices  0(t)  and  H(s)  have  been  expanded,  so  that 

* 

cancelations  for  0(t)  *  0  •  the  true  parameter  matrix,  can  occur  in  the  auxiliary 

loop  of  the  error  model  in  Fig.  4. 

From  Fig.  4,  the  error  equations  can  be  derived  as 

1 


e(t) 


ex(t) 

e2(t) 


l+cJCtKj^t) 


i+^(t)r2i£t) 


♦J(t)cct) 


♦2(tK2(t) 


where  ^(t)  -  h1;1(s)u)(t)  and  C2<t)  -  h22(s)u>(t). 
The  adaptive  laws  are  given  by 


e^t)  =  ^(t)  -  -r1e1(t)C1(t) 


e2(t)  -  ;2(t)  -  -r2e2(t)i2(t> 


I*  >0  and 


T 

r2  >  0. 


By  the  remark  following  lemma  4,  the  error  system  is  uniformly  stable  and 
e(t)  and  *(t)  are  uniformly  bounded  regardless  of  u(t) . 


*  Similar  to  the  operator  L(>)  in  [1]  for  SISO  systems,  a  diagonal  matrix  of 
operators  S(‘)  exists  which  makes  H(8)’S(s)  SPR.  For  simplicity  of  analysis 
S  (s)  ■  H(s)  is  chosen. 


(ii)  Unknown  Kp 

When  is  unknown  the  controller  structure  Is  considerably  more  involved. 
More  prior  Information  about  the  plant  transfer  matrix  is  needed  to  obtain  a 
stable  adaptive  controller. 

a)  H(s)  SPR 

The  controller  structure  remains  the  seme  as  in  Fig.  1.  The  error  model 
also  remains  the  same  as  In  Fig.  2.  The  parameter  error  matrix  *(t)  Is  updated 
according  to  the  law 

*(t)  -  6(t)  -  -re1(t)uT(t) 

The  boundedness  of  (t)  and  *(  t)  can  be  proved  by  lemma  4.  By  the  same  argument 

as  in  the  previous  case,  e^(t)  and  *(t)  -►  0  as  t  ■+■  «.  However,  now  the  input 

T 

v(t)  in  eqn.  7  Is  Kpl*(t)o>(t)  -  w  (t)r^u(t)e^(t)  ] .  It  can  be  easily  shown  that 

Kp  must  be  positive  definite  to  prove  that  V  >  0  and  V  5  0  in  lemma  4. 

For  the  transfer  matrices  in  class  I,  the  knowledge  of  the  sign  of  each 

entry  in  K  is  enough  to  meet  this  sign  definiteness  condition. 

P 

b)  H(s)  not  SPR 

Since  H(s)  is  not  SPR,  auxiliary  signals  have  to  be  generated  to  augment 

the  model  output.  To  account  for  the  unknown  gain  matrix,  an  additional  gain 

parameter  matrix  T(t)  is  introduced  in  series  with  the  auxiliary  signals.  The 

controller  structure  in  fig.  3  is  modified  slightly  to  include  ¥(t)  between  the 

summing  junctions  4  and  3.  The  error  equation  changes  somewhat  and  the  error 

model  is  shown  in  fig.  5  for  class  III  (K  nonsingular)  transfer  matrices. 

P 

Other  cases  in  classes  I  and  II  can  also  be  suitably  specialized.  ¥(0  is  also 

*  * 

adjusted  along  with  0(t)  such  that  as  t  ",f(t)  -*■  K  (K  is  a  diagonal  matrix 

* 

formed  by  the  elements  of  K  arranged  along  the  diagonal)  and  e(t)  -►  0  . 

P 

The  augmented  error  vector  from  fig.  5  for  this  case  can  be  computed  to 


Plant  Error 


V 

e2:  Auxiliary  Error 

e:  Augmented  Error 
1:  Identity  matrix  of  4v  dimension. 

Fig.  5.  Augmented  Error  Model  [unknown  Kp  and  H(s)  not  SPR] 


and 

♦^T(t)  0  “ 

c2(t) 

♦2Tco  0 

huCs)I 

53Ct) 

0  ^(t) 

h22(s)I 

V*> 

0  42T<t) 

and  G(t] 

»  dlag(g  (t)) 

-(K  )-1T(t)  -  I. 

V’1* 


h22C*)I2: 


The  corresponding  adaptive  laws  are 
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and 


®l(t)  “  'Y1ei(tKl(t)»  ®2(t)  "  _Y2ei(t)C2Ct) 

®3(t)  “  “Y3e2t3(t>’  *4(t)  *  "Y4C2(t)C4(t)* 


Yi  >  0,  1  -  1,2, 3,4. 


Using  Che  Lyapunov  function  approach  It  can  be  shown  that  the  error  system  is 
stable  and  e(t)  and  #(t)  are  bounded.  But  a  sufficient  condition  for  this  is 

that  K  be  such  that  there  exists  a  matrix  T  such  that  the  symmetric  part  of 
P 

(K  r)  be  positive  definite. 

P 

As  pointed  out  in  [2]  for  SISO  systems,  the  main  stability  question  arises 

when  H(8)  is  not  SPR.  In  this  case  auxiliary  signals  are  used  which  cannot  be 

assumed  to  be  bounded.  So  even  if  e(t)  is  bounded,  boundedness  of  the  plant 

output  y  (t)  and  hence  all  the  relevant  signals  can  not  be  concluded.  The  nature 
P 

of  the  stability  problem  is  the  same  whether  K  is  unknown  or  known.  The  fact 
.  2 

that  e(")  and  *(•)  e  L  can  be  concluded  through  the  existence  of  a  Lyapunov 

e 

function  as  was  done  in  lemma  1.  Boundedness  of  w(t)  (w(t))  end  w(t)  (u>Ct) )  can 
be  shown  using  the  same  arguments  as  in  [2]  making  it  possible  to  conclude  the 
stability  of  the  adaptive  control  loop  in  the  large. 

Remark  1 

A  similar  condition  on  the  gain  matrix  has  been  obtained  in  [8]  for  the 
discrete  case  and  in  [9]  for  the  continuous  case.  Elliott  et.al.  in  [7]  have 
made  somewhat  more  restrictive  assumptions  regarding  the  structure  of  Kp  to 
generate  stable  adaptive  laws. 
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VI.  Conclusions 

— Hermlte  normal  forms  of  nonsingular  transfer  matrices  play  a  central  role 
In  determining  the  class  of  model  transfer  matrices  which  the  plant  can  follow. 
However,  due  to  inherent  complexity  in  specifying  the  Hermlte  form  in  general, 
it  has  been  argued  that  adaptive  control  is  practically  feasible  only  for 
those  plants  which  have  diagonal  Hermlte  forms,  i.e.,  which  can  be  decoupled  by 
using  state  feedback  only.  The  sign  definiteness  of  the  high  frequency  gain 
matrix  K^j  has  been  found  to  be  sufficient  to  generate  stable  adaptive  laws. 

For  2x2  systems,  the  knowledge  of  relative  degree  of  each  scalar  transfer 
function  has  been  used  in  determining  the  Hermlte  form.  A  globally  stable 
adaptive  controller  has  been  developed  and  it  has  been  shown  that  all  2x2  stably 
invertible  systems  can  be  generically  adaptively  controlled  subject  to  the 
definiteness  condition  on  the  gain  matrix.  ^ 
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